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A rigorous QED evaluation of the two-photon exchange corrections to the g factor of lithiumlike ions is 
presented. The screened self-energy corrections are calculated for the intermediate-^ region and its accuracy 
for the high-Z region is essentially improved in comparison with that of previous calculations. As a result, 
the theoretical accuracy of the g factor of lithiumlike ions is significantly increased. The theoretical prediction 
obtained for the g factor of gth = 2.000 889 892(8) is in an excellent agreement with the corresponding 

experimental value gexp = 2.000 889 889 9(21) [A. Wagner et al., Phys. Rev. Lett. 110, 033003 (2013)]. 

PACS numbers: 31.30.J-, 31.30.js, 31.15.ac 


Highly charged ions provide not only a unique scenario for 
probing QED effects in the strongest electromagnetic fields 
but also give access to an accurate determination of funda¬ 
mental physical constants and nuclear parameters. In recent 
years, amazing progress has been made in the experimental 
and theoretical investigations of the bound-electron g factor. 
High-precision measurements of the ground state g factor of 
H-like carbon (11] and oxygen (3] and the related theoretical 
calculations provided determination of the electron mass. Re¬ 
cently, due to the substantial progress in the experimental ac¬ 
curacy of the g factor of H-like carbon and silicon the mass of 
the electron is once again substantially increased |[3. So far 
H- and Li-like silicon ions represent the heaviest ions, where 
the g factor has been measured iiS]. To date, these experi¬ 
ments provide the most stringent tests of the bound-state QED 
(BS-QED) corrections in the presence of a magnetic field. Ac¬ 
curate measurements of the g factor in few-electron ions, such 
as Li-like calcium and B-like argon (3, are already antici¬ 
pated. The investigations of the few-electron ions unlike H- 
like ions provide also an access to the many-electron QED 
corrections, which are represented by a different facet of the 
QED diagrams. 

The theoretical contributions to the g factor of Li-like ions 
can be separated into one-electron and many-electron parts. 
The one-electron terms are similar to the corresponding cor¬ 
rections to the g factor of H-like ions. The many-electron 
contributions, which define the main difference between the 
g factors of H- and Li-ions, were investigated in Refs. liQ. 
The many-electron contributions are mainly determined by 
the screened radiative and the interelectronic-interaction cor¬ 
rections. Eor low-Z ions, the screened radiative corrections 
were obtained employing the perturbation theory to the lead¬ 
ing orders in aZ (qIlIOI]- For intermediate-Z ions, the screen¬ 
ing effect was evaluated by introducing the effective screen¬ 
ing potential in the QED calculations to all orders in aZ (a. 
Eor high-Z ions, the most accurate results for the screened 
radiative corrections were obtained rigorously within a sys¬ 
tematic QED approach lEQ. The one-photon exchange 


diagrams, which represent the interelectronic-interaction cor¬ 
rection of the first order in 1/Z, were evaluated in the frame¬ 
work of QED in Ref. (3 • The second and higher-orders con¬ 
tributions of the interelectronic interaction were calculated 
by means of the large scale configuration-interaction Dirac- 
Eock-Sturm (CI-DES) method in Ref. ifioll . However, until 
now, for all values of Z the theoretical uncertainty was de¬ 
termined by the interelectronic-interaction corrections and for 
the intermediate-Z region also by the screened self-energy 
corrections. In the present Letter we report on the complete 
evaluation of the two-photon exchange and the screened self¬ 
energy corrections in the framework of a rigorous QED ap¬ 
proach within an extended Eurry picture. 


In the extended Eurry picture, to zeroth order we solve 
the Dirac equation with an effective spherically symmetric 
potential treating the interaction with the external Coulomb 
potential of the nucleus and the local screening potential ex¬ 
act to all orders. This approach significantly accelerates the 
convergence of the perturbation expansion. We use different 
types of the screening potential. The simplest choice is the 
core-Hartree (CH) potential, which is created by the charge 
density distribution of the two core electrons in the Is state. 
Other choices are the Xa potentials: Kohn-Sham (KS), Dirac- 
Hartree (DH), and Dirac-Slater (DS), which were success- 
fully employed in preyious calculations of highly charged ions 
I U i 3- 191 • Moreoyer, we haye also employed the Perdew- 
Zunger (PZ) potential lEol] and the local Dirac-Eock (LDE) 
pote ntial deriyed by inyersion of the radial Dirac equation 

M- 


Let us now turn to the eyaluation of the two-photon ex¬ 
change corrections to the g factor of Li-like ions. These cor¬ 
rections are defined by diagrams of third order in the QED per¬ 
turbation theory. The corresponding diagrams are presented 
in Eig. [T] The electron propagators in the figure haye to be 
treated in the effectiye potential (we indicate this diagrammat- 
ically yia the triple electron line). In contrast to the case of the 
original Eurry picture, in the extended Eurry picture the addi¬ 
tional counterterm diagrams appear. These diagrams are de- 
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FIG. 1: Feynman diagrams representing the second-order interelectronic-interaction corrections to the g factor in local effective potentials. 
The wavy line indicates the photon propagator and the triple lines describe the electron propagators in the effective potential. The dashed line 
terminated with the triangle denotes the interaction with the magnetic field. The counterterm diagrams are depicted in the second line. The 
symbol ® represents the extra interaction term associated with the screening potential counterterm. 


picted in the second line in Fig.[T] They are associated with an 
extra interaction term represented graphically by the symbol 
(g). Taking into account all possible permutations of the one- 
electron states, in total, we have to evaluate 36 three-electron, 
36 two-electron, and 2 one-electron diagrams, respectively. 
All together these diagrams form the complete gauge invari¬ 
ant set of the two-photon exchange contributions. The most 
difficult ones are the 16 two-electron diagrams depicted in the 
hrst line in Fig. [T] Each of this diagram contains a three-fold 
summation over the complete Dirac spectrum and an integra¬ 
tion over the loop energy. Formal expressions for the diagrams 
in the hrst line are similar to those derived for the correspond¬ 
ing calculation of the hyperhne splitting and can be found in 
Ref. II 22 I] . The formulas derived there can be taken over but, 
instead of the hyperhne-interaction potential, we employ here 
the interaction with a constant magnetic held and keep in mind 
that the Dirac spectrum is now generated by solving the Dirac 
equation with the effective potential. The derivation of the for¬ 
mal expressions for the diagrams of the second line is straight 
forward and will be presented elsewhere. More details about 
the scheme of the numerical implementations can be found 
in Ref. ll22ll . However, unlike the hyperhne splitting, in the 
case of g factor the calculations are more involved due to the 
large cancellations of various terms and poor convergence of 
the partial-wave expansion. Nevertheless, we have substan¬ 
tially increased the accuracy of the all numerical integrations 
and extended the partial-wave summation up to Kmax = 15. 
For a consistency check, we performed calculations both in 
Feynman and Coulomb gauges, and the results are found to 
be gauge invariant with a very high accuracy. 


In Table U the interelectronic-interaction corrections to the 
g factor of Li-like silicon are given. The results are ob¬ 
tained with four different starting potentials: Coulomb, core- 
Hartree, Perdew-Zunger, and local Dirac-Fock potentials. In 
the extended Furry picture, the interelectronic interaction con¬ 
tributes already in the zeroth order, due to the presence of 
the screening potential in the Dirac equation. The one-photon 


TABLE I: Interelectronic-interaction corrections to the ground-state 
g factor of Li-like ion in various starting potentials in units 

10 "®. 



Coulomb 

CH 

PZ 

LDF 

Zeroth order 


348.267 

321.632 

349.636 

First order 

321.592 

-33.549 

-5.990 

-33.846 

Second order 

-6.876 

0.137 

-0.866 

-0.976 

Higher orders 

0.085(22) 

-0.046(6) 

0.034(6) 

-0.005(6) 

Total 

314.801(22) 

314.809(6) 314.810(6) 

314.808(6) 


(hrst order) and two-photon (second order) exchange correc¬ 
tions have been evaluated to all orders in aZ in the frame¬ 
work of rigorous QED approach. The higher-order correc¬ 
tions have been extracted from the calculations performed 
by means of the large-scale conhguration-interaction Dirac- 
Eock-Sturm method described in Refs. M As it was 
expected, the employment of the extended Eurry picture in¬ 
creases the convergence of the perturbation expansion. This 
allows us to reduce the absolute uncertainty of the higher or¬ 
der interelectronic-interaction corrections. Einally, the rigor¬ 
ous evaluation of the two-photon exchange corrections and the 
improved calculations of the higher-order terms allow us to 
signihcantly increase the total accuracy of the interelectronic- 
interaction terms for all ions under consideration. E.g., in the 
case of ion the previous result was 0.000 314 903(74) 

ifl^ . while the present calculation yields 0.000 314 809(6), 
and in the case of 208p|^79-i- instead of previous value 
0.002 140 7(27) ifl^ we now receive 0.002 139 34(4). 

Let us now turn to the screened self-energy corrections to 
the g factor of Li-like ions. In Refs. lll2l 1 1311 these correc¬ 
tions have been rigorously evaluated only for the 208p|^79-i- 
and 238pj89-i- ions. The reasons for this are twofold. Eirst is 
the large numerical cancellations which occur in the point- 
by-point difference. Second is the poor convergence of the 
partial-wave expansion. In order to overcome these problems, 
we have performed the calculations in the extended Furry pic- 















































3 


ture and employed a special treatment of the many-potential 
terms. The Feynman diagrams in the extended Furry pic¬ 
ture corresponding to the screened self-energy corrections to 
the g factor are presented in Fig. jH The corresponding ex¬ 
pressions derived in Refs. il2lll3n remain formally the same 
but keeping in mind that the Dirac spectrum is now gener¬ 
ated by solving the Dirac equation with the effective poten¬ 
tial. In the second line of Fig. |2] the additional counterterm 
diagrams are depicted. The derivation of the formal expres¬ 
sions for them is relatively simple and will be presented else¬ 
where. The employment of the extended Furry picture al¬ 
lows us to substantially reduce the numerical cancellations 
of different terms as well as to improve convergence of the 
partial-wave expansion. However, in order to improve the 
convergence even further, we have employed a specific treat¬ 
ment of some many-potential terms. The standard way to 
treat the vertex and reducible corrections is to separate terms 
(zero-potential contributions) in which bound-electron propa¬ 
gators are replaced by free propagators. The remaining many- 
potential terms being ultraviolet finite are generally calculated 
directly in coordinate space ll24ll . However, for gaining bet¬ 
ter control over the partial-wave summation we separate also 
the so-called one-potential contributions. In this way the one- 
potential terms are treated in the momentum space. Such 
treatment of the one-potential term was applied in previous 
calculations of the one-electron self-energy corrections to the 
29|] and to the magnetic-dipole transi- 
l30l] . Here, we extend this procedure 


g factor in Refs. |l25 
tion amplitude in Ref. 
to the evaluation of the screened self-energy corrections to the 
g factor. Performing the analysis of the convergence of the 
partial-wave expansion for different terms we have found, that 
such treatment should be applied to the terms (Cl) Eq. (32), 
(H3) Eq. (38), (II) Eq. (47), and (13) Eq. (49) in Ref. |[I1. The 
corresponding one-potential contributions are given by the ex¬ 
pressions: 
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xV'r,„(p') + (« ^ ^): (2) 


where p = {Sa, p), p' = {Sa, p')’ d = {Sa, q), A=£a- £pa, 
and the notation (a O b) stands for the contribution with in¬ 
terchanged labels a and 6; 7 ^ = ( 70 , 7 ) are the Dirac matri¬ 
ces, S'f(p) = (7 • P — m)~^ is the free-electron propagator, 
the interelectronic-interaction operator /(e) and its derivatives 
are defined in a similar way as in Ref. joll . and 14ff is the ef¬ 
fective potential being the sum of the nuclear and screening 
potentials. Tq is the operator of interaction with a constant 
magnetic field, which reads in the momentum space: 

To(p) = i/ro(27r)3[a x Vp(5^(p)] • "H , (5) 

where po = | e| /2 is the Bohr magneton and H is the magnetic 
field directed along the z axis. The wave function \pa) is given 
by the expression: 

\Va) = {Cb\PaPb\To\a) + {Ca\PbPa\To\b) 

+ {ab\r{A)\PaPb)(^{a\To\a) - ^{b\To\b)) 

+ \^a){ci-b\I{A)\PaPb) + |Cb|PaPh)(o|To|o) I J (6) 

and the wave functions |^) and |(^) are defined similar as in 
Ref. d. 

The ultraviolet-finite one-potential terms given by Eqs. (1)- 
(4) have been evaluated in the momentum space. The corre¬ 
sponding expressions in the coordinate space have been sub¬ 
tracted from the related many-potential terms by means of 
point-by-point difference. The partial-wave expansion for the 
many-potential terms was terminated at Kmax = 15, and the 
remainder of the sum was estimated by a least-square polyno¬ 
mial fitting and by the e-algorithm of the Fade approximation. 
As a result, we have significantly increased the accuracy of 
the screened self-energy correction. In the case of 2®Si^^+ ion 
the previous result was —0.000 000 218(46) ifl^ . while the 
present calculation yields —0.000 000 242(5), and in the case 
of 208p|^79+ jQp instead of previous value —0.000 003 3(2) 
idllH] we now receive —0.000 003 44(2). 

In Table im the individual contributions and the total values 
of the g factor for Li-like silicon calcium ^°Ca^^+, 






















4 



FIG. 2: Feynman diagrams representing the screened self-energy corrections to the g factor in local effective potentials. The wavy line 
indicates the photon propagator and the triple lines describe the electron propagators in the effective potential. The dashed line terminated 
with the triangle denotes the interaction with the magnetic field. The counterterm diagrams are depicted in the second line. The symbol (g) 
represents the extra interaction term associated with the screening potential counterterm. 


lead 208p^79-i- uranium 238pj89+ presented together 
with the previously reported theoretical results and the exper¬ 
imental value for the case of silicon. The screened self-energy 
and interelectronic-interaction corrections calculated in this 
Letter allow to substantially increase the theoretical accuracy 
for all ions under consideration. The other contributions to the 
g factor presented in Table HJ were considered in detail in our 
previous studies CMl- Comparison with the experimental 
value for Li-like silicon ion provides tests of relativistic in- 
terelectronic interaction on a level of 10“®, the one-electron 
BS-QED on a level of 0.7%, and the screened BS-QED on a 
level of 3%. Thus, the current studies provide the most ac¬ 
curate test of the many-electron QED effects in the case of g 
factor. The further improvement of the g factor theory for Li- 
like ions requires at first the rigorous evaluation of the three- 
photon exchange diagrams and the subsequent betterment of 
the screened self-energy contribution for the intermediate-Z 
region, and the one-electron two-loop and nuclear recoil cor¬ 
rections for the high-Z region. 

The techniques and numerical methods developed can also 
be extended for the g factor of B-like ions, where the corre¬ 
sponding studies can also lead to an independent determina¬ 
tion of the fine-structure constant ini]. 
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TABLE II: Individual contributions to the ground-state g factor of Li-like ions and comparison with the previously reported theoretical values 
as well as with the experimental result for the ion. 



g8gjii+ 


2U8p^jV9-l- 

1^38-^89+ 

Dirac value (point nucleus) 

1.998 254 751 

1.996426 011 

1.932 002 904 

1.910722 624 

Finite nuclear size 

0.000000 003 

0.000 000014 

0.000078 57(14) 

0.00024162(36) 

QED, ~ a 

0.002 324 044 

0.002 325 555(5) 

0.0024117(1) 

0.0024463(2) 

QED, ~ 

-0.000003 517(1) 

-0.000 003 520(2) 

-0.000003 6(5) 

-0.000003 6(8) 

Interelectronic interaction 

0.000 314 809(6) 

0.000454290(9) 

0.002139 34(4) 

0.002 500 05(6) 

Screened self-energy 

-0.000000 242(5) 

-0.000 000387(7) 

-0.000003 44(2) 

-0.000004 73(3) 

Screened vacuum-polarization 
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0.000001 53(3) 
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Nuclear polarization 



-0.00000004(2) 

-0.000000 27(14) 

Total theory 

2.000 889 892(8) 

1.999 202 041(13) 

1.936627 2(6) 

1.915904 8(11) 


2.000 889 909(51)“ 
2.000 890 005(87)*’ 

1.999 202 24(17)*' 

1.936628 7(28)“ 

1.915 905 7(41)“ 

Experiment 

2.000 889 889 9(21)“ 
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